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Abstract 

We consider a generalization of the Lotka-McKendrick problem describ- 
ing the dynamics of an age-structured population with time-dependent vi- 
tal rates. The generalization consists in allowing the initial and the bound- 
ary conditions to be derivatives of the Dirac measure. We construct a 
unique D'-solution in the framework of intrinsic multiplication of distribu- 
tions. We also investigate the regularity of this solution. 
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1 Introduction 

We consider a non-classical hyperbolic problem with integral boundary condition 

{d t + d x )u = p(x,t)u + g(x,t), (x,t)eU (1) 
u\ t=0 = a(x), x e [0,L) (2) 

L 

u\ x= o = c{t) J b{x)udx, te [0, oo) , (3) 



where 



n = {(x, t) G IT | < x < L, t > 0}. 



From the point of view of applications, (l)-(3) describes the dynamics of the age- 
structured population (see i.e. [1, 3, 15, 23, 28]). There u denotes the distribution 
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of individuals having age x > at time t > 0, a(x) is the initial distribution, 
—p(x, t) denotes the mortality rate, b(x) denotes the age-dependent fertility rate, 
c(t) is the specific fertility rate of females, g(x,t) is the distribution of migrants, 
L is the maximum age attained by individuals. Furthermore, b(x) = on [0, L\ \ 
[Li,L 2 ], where [L 1: L 2 ] C [0, L\ is the fertility period of females. The evolution 
of u without diffusion is governed by (l)-(3). The system (l)-(3) is a continuous 
model of a discrete structure. As in many problems of such a kind, it is natural 
to consider singular initial and boundary data. We focus on the case when these 
data have singular support in finitely many points, i.e. 



a(x) = a r (x) + ^ du^ mi \x — Xi) for some du £ R, TOj £ N , x { £ (0, L), 

s 

b(x) = b r (x) + d 2i 5 {nk \x - x k ) for some d 2i £ M, n k £ N , x k £ (0, L), (4) 
fc=i 

q 

c(t)=c r (t)+Y J d^\t-t j ) for some d 3i £ R, lj £ N , tj £ (0, oo). 

3=1 

The data of the Dirac measure type enable us to model the point-concentration 
of various demographic parameters. 

The problem under consideration is of interest from both biological and math- 
ematical points of view. 

Biological motivation A basic model describing the evolution of an age- 
structured population is given by the Lotka-McKendrick system: 

(d t + d x )u = -p(x)u (5) 
u\ t=0 = u (x) (6) 



u\ 



= o = J b(x)udx. (7) 



The differential equation describes the aging of the population and the output 
due to deaths. The integral u(x, t) dx gives the number of individuals at time 
t having age x in the range d\ <J x ^ ct 2 . Thus, the third equation is responsible 
for newborns, entering the population at age zero. 

A biological generalization of (5) to (l)-(3) consists in the allowing the fertility 
and mortality rates to depend on t (see e.g. [9, 10, 14]). In reality the vital rates 
are never time-homogeneous and adapt to the changing social and technological 
environment. Introducing the (^-distributional data in (2) and (3) also has a 
biological meaning (see [15]). 
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In demography, c(t) is the total fertility rate of the population at time t, in 
other words, the average number of childbirths per female during her reproductive 
period. On one side, the results presented in the paper could shed a new light 
on the so-called c-control problems when one wants to control the population 
only through changing c(t). Chinese scientists used discrete models to provide 
mathematical background for the unicity child policy (c-control problem) in the 
People's Republic of China [25, 26, 29]. Continuous models in the context of the 
c-control problem were considered in [8]. In contrast to the aforementioned pa- 
pers, the presence of strongly singular data in (2) and (3) allows one to combine 
the continuity of the model with the discreteness of the real evolutionary process. 
Occurrence of strong singularities in c(x) can be motivated by synchronized and 
concentrated reproduction of the species. This also allows one to involve statisti- 
cal data into (l)-(3) and perhaps makes our model competitive with discrete-time 
and discrete-age models [2]. 

Involving strong singularities into the model could have another interpreta- 
tion: such singularities can be produced by a linearization of nonlinear problems 
with discontinuous data. Thus this opens a space for interesting nonlinear con- 
sequences. 

Mathematical motivation We consider our paper as a further step in the 
study of generalized solutions to initial-boundary hyperbolic problems in two 
variables. 

Since the singularities given on dU expand inside II along characteristic curves 
of the equation (1), a solution preserves at least the same order of regularity as it 
has on 911. This causes multiplication of distributions under the integral sign in 
(3). In spite of this complication, we find distributional solutions of (l)-(3). In 
parallel, we study propagation, interaction and creation of new singularities for 
the problem (l)-(3). 

Initial-boundary semilinear hyperbolic problems with distributional data were 
studied, among others, in [18, 11, 12]. There also appears a complication with 
multiplication of distributions that is caused by nonlinear right-hand sides of 
the differential equations and also by boundary conditions that are nonlinear 
(with bounded nonlinearity) in [18], nonseparable in [12], and integral in [11]. To 
overcome this complication, the authors use the framework of delta waves (see 
[20]). In other words, they find solutions by regularizing all singular data, solving 
the regularized system and then passing the obtained sequential solution to a 
weak limit. 

Boundary and initial-boundary value problems for a linear second order hy- 
perbolic equation [22] and the general strictly hyperbolic systems in the Leray- 
Volevich sense [21] are studied in a complete scale of Sobolev type spaces depending 
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on parameters s and r, where s characterizes the smoothness of a solution in all 
variables and r characterizes additional smoothness in the tangential variables. 
Sobolev-type a priori estimates are obtained and, based on them, the existence 
and uniqueness results in Sobolev spaces are proved. 

In contrast to the aforementioned papers we here treat integral boundary 
conditions and show that the problem (l)-(3) is solvable in the distributional 
sense. We construct a unique distributional solution by means of multiplication 
of distributions in the sense of Hormander [7] . 

We show that the boundary condition (3) causes anomalous singularities at 
the time when singular characteristics and vertical singular lines arising from the 
data of (3) intersect. In the case that the singular part of b(x) is a sum of deriva- 
tives of the Dirac measure, the solution becomes more singular. In the case that 
the initial and the boundary data are Dirac measures, the solution preserves the 
same order of regularity. Similar phenomenon was shown in [27] for a semilinear 
hyperbolic Cauchy problem with strongly singular initial data, where interaction 
of singularities was caused by the nonlinearity of the equations. Anomalous sin- 
gularities were considered also in [19] and [17], where propagation of singularities 
for, respectively, initial and initial-boundary semilinear hyperbolic problems were 
studied. There was proved that, if the initial data are, at worst, jump disconti- 
nuities, then the singularities at the common point of singular characteristics of 
the differential equations are weaker. Furthermore, if boundary data are regular 
enough, then reflected singularities cannot be stronger than the corresponding 
incoming singularities. It turns out [4, 13] that in some cases of nonseparable 
boundary conditions the solution becomes more regular in time, namely for C 1 - 
initial data it becomes fc-times continuously differentiable for any desired k EN 
in a finite time. 

Organization of the paper Section 2 contains some basic facts from the 
theory of distributions. In Section 3 we describe our problem in detail and state 
our result. Sections 4-9 present successive steps of construction of a distributional 
solution to the problem. In particular, the integral boundary condition is treated 
in Section 5. In parallel we analyze the regularity of the solution. The uniqueness 
is proved in Section 10. 

2 Background 

For convenience of the reader we here recall the relevant material from [5, 6, 7, 24] 
without proofs. Throughout the paper we will denote by (•,•): V x V — > R the 
dual pairing on the space V of C°°-functions having compact support. 

Definition 1 ([6], 2.5) A distribution u G V(R 2 ) is microlocally smooth at 
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(x, t,£,r)) ((C,v) 0) if the following condition holds: If u is localized about (x,t) 
by ip G V(R 2 ) with tp = 1 in a neighborhood of (x,t), then the Fourier transform 
of tpu is rapidly decreasing in an open cone about (£,77). The wave front set of u, 
WF(u), is the complement in M 4 of the set of microlocally smooth points. 

Proposition 2 ([7], 8.1.5) Let u G "D'(lR 2 ) and P(x,D) be a linear differential 
operator with smooth coefficients. Then 

WF(Pw) C WF(m). 

Definition 3 ([7], 6.1.2) Let X, Y C K 2 be open sets and u G V'{Y). Let f : 
X — > Y be a smooth invertible map such that its derivative is surjective. Then 
the pullback of u by /, f*u, is a unique continuous linear map: V'{Y) — > V{X) 
such that for all ip G T>(Y) 

(r»M = (u,w- 1 )\{<pof- 1 )), 

where is the Jacobian matrix of 

Theorem 4 ([7], 8.2.7) Let X be a manifold and Y a submanifold with normal 
bundle denoted by N(Y). For every distribution u in X with WF(-u) disjoint from 
N(Y), the restriction u\y of u to Y is a well-defined distribution on Y that is the 
pullback by the inclusion Y <^-> X. 

Theorem 5 ([7], 5.1.1) For any distributions u G V'(Xi) and v G T>'(X 2 ) there 
exists a unique distribution w G T>'(Xi x X 2 ) such that 

(w, (pi <S> ¥2) = (u, (pi) (v, ip 2 ) , (pi G T>(Xi) 

and 

(w, (p) = (u, (v, (p(x 1 ,x 2 ))) = (v, (u, (p(xi, x 2 ))) , ip G V(X 1 x X 2 ). 

Here u acts on ip(xi,x 2 ) as on a function of Xi and v acts on ip(x 1 ,x 2 ) as on a 
function of x 2 . 

The distribution w as in Theorem 5 is called the tensor product of u and v, and 
denoted by w = u <g> v. 

Theorem 6 ([5], 11.2.2) Let X,Y be open sets in M 2 and let f : X — > Y be a 
diffeomorphism. If u G V'(Y), then f*u, the pull-back of u, is well defined, and 
we have 

WF(/») = {(x,df x r]) : (f(x),r)) G WF(«)}. 



5 



Theorem 7 ([7], 8.2.10) Ifv,w G V'(X), then the product v ■ w is well defined 
as the pullback of the tensor product v <E> w by the diagonal map 5 : R — > R x R 
unless (x, t, £, rf) G WF(t>) and (x, t, — £, —77) G WF(iu) /or some (x, t, £, rj) . 

Theorem 8 ([24], 8.6) If a distribution u is identically equal to on each of the 

domains G i} i > 1, then u is identically equal to on G = U G,. 

i>i 

3 Statement of the results 

For simplicity of technicalities we assume that both the initial and the boundary 
data have singular supports at a single point and are the Dirac measures or 
derivatives of the Dirac measure. This causes no loss of generality for the problem 
if the singular parts of the initial and the boundary data are finite sums of the 
Dirac measures and derivatives thereof, i.e. they are of the form (4). Specifically, 
we consider the following system 

p(x,t)u + g(x,t), (8) 
a r (x) + 5 {m \x-x* 1 ), xe[0,L) (9) 

L 

fait) + 6<i\t - tj) J{b r {x) +5^{x-x 1 ))udx, (10) 



t e [0, 00), 

where x\ > 0,xl > 0,£i > 0, and m,j,n G N . Without loss of generality we can 
assume that x\ < x\. 

We impose the following conditions: 

Assumption 1. a®(0) = 0, c^(0) = for all i G N . 

Assumption 2. bfi(L) = for all i e No and there exists e > such that 
b r (x) = for x G [0, e\. 

Assumption 3. The functions p and g are smooth in R 2 , a r is smooth on [0, L), 
b r is smooth on [0, L], and c r is smooth on [0, 00). 

Note that Assumption 1 ensures the arbitrary order compatibility between 
(9) and (10). Assumption 2 not particularly restrictive from the practical point 
of view, since [0, L\ covers the fertility period of females. 

All characteristics of the differential equation (1) as solutions to the following 
initial problem for ordinary differential equation: 

dx 

— = 1, x(t ) = x , 
where (xq, to) G R 2 , are given by the formula x = t + x$ — t . 



(d t + d x )u = 
u\ t =o = 

u\ x=0 = 
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Definition 9 Let I + = \J n>0 I + [n], where I+[n] are subsets of R 2 defined by in- 
duction as follows. 

• I + [0] is the union of the characteristics x = t + x\ and x = t — t\ . 

• Let n > 1 . If I + [n — 1] includes the characteristic x = t + X\ — i , then I + [n] 
includes the characteristic x — t — t. 

For characteristics contributing into /+ denote their intersection points with the 
positive semiaxis x = by i*,^ ■ ■ • • We assume that t* < t* +1 for j > 1. The 
union of all singular characteristics of the initial problem, as it will be shown, is 
included into the set I + . In fact, we will show that singsuppw C i+. 

Assumption 4- %\ — t\ ^ x\, t x — x\ ^ t* for all t* < t±. 

This assumption excludes the situation when three different singularities in- 
tersect at the same point. Without this assumption the distributional solution 
does not exist, because there appears multiplication of the Dirac measure onto 
itself. 

Our goal is, using distributional multiplication, to obtain distributional solu- 
tion to (8)-(10). We use the notion of the so-called " WF favorable" product which 
is due to L. Hormander [7] and is in the second level of M. Oberguggenberger's 
hierarchy of intrinsic distributional products [16, p. 69]. 

We actually obtain distributional solution in a domain Q C M 2 that is the 
domain of influence (or determinacy) of the problem (8)-(10). Clearly, Q is the 
union of all characteristics x = t + x — t passing through those points (x ,t ) 
on the boundary of II where the conditions (9) and (10) are given, i.e. through 
points (x ,t ) e ([0, L) x {0}) U ({0} x [0, oo)). In other words, 

Q = {(x,t) G R 2 \x < t + L}. 

Definition 10 A distribution u is called a V'{VL)- solution to the problem (8)- 
(10) if the following conditions are met. 

1. The equation (8) is satisfied in V'(Q): for every <p e V(Q) 

((d t + d x - p(x, t))u, ip) = (g(x, t),<p). 

2. u is restrictable to [0, L) x {0} in the sense of Hormander (see Theorem 4) 
and u\ t =o = a r (x) + 5^ m \x — x\), x G [0, L). 

3. The product of (b r (x) + 5^ n \x - x 1 )) <S> l(t) and u(x,t) exists in T>'(U) in the 
sense of Hormander (see Theorem 7). 

4- Jo (br(x) + 5( n \x — xifj ® l(t) udx is a distribution v G V'(R + ) defined by 
(v,ij(t)) = ([(b r (x) + 5 {n) (x-x 1 ))®l(t)]u,l(x)^ij(t)), rj){t) £%), 
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where b r (x) = 0, x [0,L]. 

5. v is a smooth function in t\ . 

6. u is restrictable to {0} x [0, oo) in the sense of Hormander (see Theorem 4) 
and u\ x=0 = (c r (t) + 5^\t — ti))v, t E [0, oo). 

7. sing supp u C \ {{x, t) \ x — t} . 

Our next objective is to define the solution concept for (8)-(10) on II. It is not 
so obvious how we should define the restriction of u e V(Tl) to the boundary of 
IT so that the initial and the boundary conditions are meaningful. In this respect 
let us make the following observation. 

Note that TT\{(L,0)} C Q. Let Q C D, be a domain such that TT\ {(L, 0)} C 
Qo and u be a X > '(Jl)-solution to the problem (8)-(10) in the sense of Definition 10. 
Then u restricted to VLq is a ©'(f^-solution to the problem (8)-(10) in the sense 
of the same definition. This suggests the following definition. 

Definition 11 Let u be a V (Vl) -solution to the problem (8) -(10) in the sense 
of Definition 10. Then u restricted to II is called a V '(II) -solution to the prob- 



We are now prepared to state the existence result. 

Theorem 12 1. Let Assumptions 1-4 hold. Then there exists a V{VL) -solution u 
to the problem (8) -(10) in the sense of Definition 10 such that: 



2. Let Assumptions 1~4 hold. Then there exists a V (II) -solution to the problem 
(8) -(10) in the sense of Definition 11. 

Given a domain G, set 



lem (8)-(10). 



Set 



n + = {( x ,t) e n\x > o,t > o}. 



the restriction of u to any domain Q! + D fi+ such that any 
characteristic of (8) intersects dQ' + at a single point does not 
depend on the values of the functions p and g on Q\ Q', . 



(11) 



V' + {G) = {ue V(G) 



u = for (x < or t < 0)} . 



Write 



FK ' 1 [o if (x,t) e n\n+. 

Similarly to p we define a modification of g and denote it by g. 
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Definition 13 u G V + (Q) is called a V' + (Vt)- solution to the problem (8)-(10) if 
the following conditions are met. 

1. Items 3-5 of Definition 10 hold. 

2. The equation (8) is satisfied in V + {VL): for every <p G V{VL) 



d t + X(x, t)d x - p{x, t))u, if) = (g(x, t), (p) 
+ ((a r (x) + 5 (m \x - x\)) <g> S(t) + 5(x) <g> [(c r (t) + 6®(t - h))v] : <p) , 



where a r (x) =0 if x < and v(t) = ift<0. 
3. sing supp u \ dVL + C Q + \ { (x, i) \ x = t } . 

Proposition 14 Let u be a V (VL)- solution to the problem (8)-(10) in the sense 
of Definition 10 that satisfies (11). Then there exists a V + (Q)- solution u to the 
problem (8)-(10) in the sense of Definition 13 such that 

u = u in V'(Q + ). 

This proposition is a straightforward consequence of Definitions 10 and 13. Since 
II C Q + , it makes sense to state the uniqueness result in V' + (fl). Write 



S(x, t) = exp < 



J p(r;x,t,T)dT 



\0(x,t) 



(12) 



where 9(x,t) — (t — x)H(t — x), H(z) is the Heaviside function. We write S for 
the function S given by (12), where p is replaced by —p. 

Assumption 5. For every T > there exists T > T such that S(x, T) ^ 
for all x such that (x, T) G Vt + . 

Theorem 15 1. Let Assumptions 1-5 hold. Then a V' + (Q) -solution to the prob- 
lem (8)-(10) is unique. 

2. Let Assumptions 1-5 hold. Then a V (II) -solution to the problem (8)-(10) is 
unique. 

From the construction of a ©'(f^-solution presented in the proof of Theo- 
rem 12 we will see that in general there appear new singularities stronger than 
the initial singularities. In other words, the singular order (cf. [24, §13]) of the 
distributional solution grows in time. We state this result in the following theo- 
rem. 
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Theorem 16 1. Let u be the V(U) -solution to the problem (8)-(10), where 
n > 1 and S(x,t) ^ for all (x,t) G IT. Then for each % > 1 there exist j > i 
and n' > 1 suc/i t/iat the singular order of u is equal to n' in a neighborhood of 
x — t — t* and the singular order of u is equal to n' + n in a neighborhood of 
x = t-t*. 

2. If n — j — m — 0, then the singular order of u on II is equal to 1. 

We now start with the proof of Theorem 12 which will take Sections 4-9. It 
is sufficient to solve the problem in the domain 

Q T = {( x , t) e D, 1 1 - T < x, -T < t < T} 

(see the picture) for an arbitrary fixed T > 0. Observe that f2 T is the intersection 
of the strip E x (— T, T) with the domain of determinacy of (8) with respect to 
the set ([0,L) x {0}) U ({0} x [0,T)). Fix T > 0. We start with a subdomain 

n% = {( x , t) e tt T 1 1 < x < t + L}. 

4 The solution on 

Observe that fig is the intersection of the strip R x (— T, T) with the domain of 
determinacy of the problem (8)-(9). In the case that the initial data are functions, 
a unique solution to the problem (8)-(9) on f2g can be written in the form 

u(x, t) = Si(x, t) + S(x, t)a r (x — t) + S(x,t)5 {m) (x-t-xl) (13) 

with the functions S(x,t) given by (12) and 

t 

J p(T + x — t),Tdr 

J piji + x — t, Ti) dri > g{r + x — t, t) dr. (14) 

x,t) J 

Let Ai(x,t) = 5^(x) ® l(t) and Bi(x,t) = l(x) (g>5 (i) (t) be the distributions in M 2 
that are derivatives of the Dirac measure S^(x) and S^(t) supported along the 
t-axis and the x-axis, respectively. They are defined by the equalities 




(A(x,t), v ( x ,t)) = (-iyJ<pW(o,t)dt, 



ii 



{Bi(x,t),<p(x,t)) = (-iyj <p?\x,0)dx 



for all tp G V(R 2 ). When % — 0, then we have the Dirac measure supported along 
the respective axes. 

Let / be the smooth map 

/ : (x,t) -> (x,x - t - xl). 

The inverse 

f- 1 : (a:,*) -> (x,x-t-x\) 

is unique and maps the rr-axis to the curve t = x — x\ and the t-axis onto itself. 
Moreover, 

/'0M)=(1 _°!). 

Hence the Jacobian of / 

J(f) = |/'| = -1 ^ 

and /*.B m = — £ — #*), the pullback of B rn by / (see Definition 3), is 

well defined. Therefore the distribution 5^ m \x — t — x\) acts on test functions 
(p E V{R 2 ) in the following way: 



= | d?<p{x,x- t - xt)\ t=Q dx = -j dr^(x,t)\ t=x _ xt dx. 

Hence, similarly to B m , f*B m is the m-th derivative of the Dirac measure sup- 
ported along the line t = x — x\. 

Definition 17 A distribution u is called a V '(Qq) -solution to the problem (8), 
(9) if Items 1 and 2 of Definition 10 with Q replaced by Qq hold. 

Lemma 18 u(x,t) given by the formula (13) is a V {VL^)- solution to the problem 
(8)-(9). 

Proof. A straightforward verification shows that the sum of the first two 
summands in (13) is a smooth (and, therefore, distributional) solution to the 
problem (8)-(9) with the singular part of the initial condition (9) identically equal 
to 0. Our goal is now to prove that the third summand in (13) is a distributional 
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solution to the homogeneous equation (8) with singular initial condition 5^ (x — 
x*). Indeed, for all (p G V(Hq), we have 

((dt + d x )(S8( m \x-t-xl)),<p) 

= _(S5( m \x-t-xl),d t <p + d x ip) 

= - (5 {m \x -t-x\), Sd t cp + Sd x <p) 

= _ ^M( x _ t _ x *), dt ( S(p ) + d x (S<p) - d t S<p - d x Sip) . 

Since w = 5^ (x — t — x*) is a distribution in x — t, this is a weak solution to the 
equation (d t + d x )w = 0. Note that Sip G V(Qq). Therefore 

( S W(x - t - xl), d t (S<p) + d x (Sip)) = 0. 

By (12), d t S + d x S = pS. The desired assertion is therewith proved. 

It remains to prove that S{x,t)5^ m \x-t-x\) may be restricted to the initial 
interval X = [0, L) x {0}. For this purpose we use Theorems 4 and 6. Observe 
that / restricted to Qq is a diffeomorphism. We check the condition 

WF(Sf*B m )nN(X) = <b, (15) 

where the normal bundle N(X) to X is defined by the formula 

N(X) = {(x, t, £, 77) | (x, t) G X, (T {x>t) (X), (£, 77)) = 0} 

and T( Xjt )(X) is the space of all tangent vectors to X at (x,t). It is clear that in 
our case 

iV(X) = {(x,0,0,77),77 ^0}. 
Let us now look at WF(Sf*B m ). By Proposition 2, we have 

WF(Sf*B m ) C WF(f*B m ). 

Recall that by definition 

WF(/'fl m ) = {(x : t : df x • (£,77)) : (f(x,t),t,r,) G WF(B m )}. 
We also have 

WF(B m ) C WF(B ) = {(x, 0, 0, 77), 77 ^ 0}. 

It follows that f(x,t) is equal to (x, 0). Therefore (x, t) = (x, x — #*). Further- 
more, 

^=(J -l)' ^•(°^)=(_ ?7 )- 
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As a consequence, 

WF(Sf*B m ) C {(x,x-xl,r), -r)),r) + 0}. 

This means that S(x,t)5 (m) (x — t — xl) is restrictable to X. Considering the 
distribution 5^ m \x — t — x\) to be smooth in t with distributional values in x, the 
initial condition (15) follows from (13). This finishes the proof. □ 

We have proved that u defined by (13) satisfies Items 1 and 2 of Definition 10 
with Q replaced by Qq. Items 4-7 on Qq do not need any proof. Item 3 will be 
given by Lemma 20 in the next section. 



5 Multiplication of distributions under the in- 
tegral in (3) 

In the further sections we will extend the solution over 

SI* = t) e tt T 1 1 - T < x < t}. 

We use the fact that any D'(fi) -solution u to our problem is representable as 

u(x, t) — u (x, t) + Ui(x, t), (16) 

where u = u in V'(Qq), u is identically equal to on Qf , u± — u in V'(Qj), and 
Mi is identically equal to on Qq . Indeed, if u is a solution, then it is a smooth 
function in a neighborhood of {(x,t)\x = t} (see Item 7 of Definition 10). For 
an arbitrary (p G V(Q T ) consider a representation (p(x,t) = ipi(x,t) + (p 2 {x,t) + 
ip 3 (x,t) such that ipi(x,t) G V(Q T ), supp^i C Qq , supp ip 2 H sing supp u = 0, and 
suppy? 3 C Qj. Hence 

(u + Mi, ip) = (u , tfi + tf 2 ) + (til, ^2 + tf3) 

= (u, y?i) + (u , (f 2 ) + (til, f2) + (ti, ^3) = (ti, + + <f3) = (ti, p) ■ 

Using (16), we rewrite v(t) (see Item 4 of Definition 10) in the form: 

L L 

v(t) = J b(x)uo(x, t) dx + J b{x)u\{x,t) dx. 


In this section we compute the integral 

L 

I (t) = J b(x)u (x, t)dx, 0<t< T, (17) 
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that will be used in the construction. We have to tackle the multiplication of 
distributions involved in the integrand. For technical reasons we extend a r (x) 
and b r (x) over all R defining them to be outside [0, L}. By (13), we rewrite (17) 
as follows 



L 



I (t) — J b r (x)(S(x, t)a r (x — t) + Si(x, t)) dx 
t 

L 

+ J 5 {n) (x-x 1 )(S(x,t)a r (x-t) + S 1 (x,t))dx 
o 

L 

+ J b r {x)S{x,t)5 {m) {x-t-xl)dx 



L 

+ f 5 {n) {x - x 1 )S{x,t)5 {m) (x - t - xl) dx. 



To compute the second integral we take a test function ip{t) e £>(0,T) and 
consider the dual pairing (see Definition 10, Item 4) 

(5 {n \x - x 1 )(S(x,t)a r (x - t) + l(x) ® ip{t)) 

= (5 {n \x - Xl ) <g> 1(0, (S(x, t)a r (x - t) + (a, t))^(t) 
= (-1)" (l(t), (S(x,t)a r (x -t) + S 1 (x,t))^\ x=Xl ij(t)) 
= (-1)" ((S(x, t)a r (x -t) + S 1 (x, t))W \ x=xi ,m 



Let us compute the third integral: 

(S(x, t)b r (x)5 im \x - t - xl), l(x) <g> i/>(t)] 

= (q*5^(x),S(x,t)b r (x)^(t)) 
= (5 {m \x), (S(x + t + x* 1: t)b r (x + t + x\)^(t)) o ^ 

= (-ir sr (^(x + 1 + xi,t)b r ( X + 1 + xi)) \ x=o , m] 

= (-l) m (S(x + t + xl,t)b r (x + t + x{) , V(t)) • 

\ 21 = / 

To compute the last integral in the expression for I (t) we need the following fact. 

Lemma 19 The product of two distributions v = S^ n \x — X\) ® l(t) and w = 
§(rn)^ x _i _ ex ists in the sense of Hormander (see Theorem 7). 
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Proof. Recall that 

WF(.) = KUi,OUi^O}, 

WF(w) c {x,x - xl&fa-b^O}. 
Thus all conditions of Theorem 7 are true and the lemma follows. □ 
We have proved the following lemma. 

Lemma 20 A distribution u defined by (13) satisfies Item 4 of Definition 10 with 
II replaced by II n VLq . 

Turning back to computing the last integral in Io(t), consider the map 

H : (x,t) — > (x — Xi, x — t — x\) 

and the inverse map 

H^ 1 : (x, t) — > (x + xi, x — t + Xi — x\). 

Define H*A n = 5^ n \x - x x ) ® l(t) and H*B m = 5^ m \x -t-x\). Let us check 
that the former definition is unambiguous: For any tp G V(M. 2 ) we have 

(H*A n , ip(x, t)) = - (A n , ip(x + x ± , x - t + x ± - z*)) 
= ~(^5 (n) (x), J ip(x + xi,x -t + Xi -xl))dtj 

= ^ n) (x),J <p(x + x u r) dr^i 

= (-l)"| <p<?\x 1 ,T)dT=(6<- n \x-X 1 )®l(t)Mx,t)). 

Here we used a simple change of coordinates t — > t — x — t + Xi — x\. 

We are now in a position to compute the product of two distributions 5^ (x — 
x±) and 5( m \x — t — x\): For any (p G V(R 2 ) we have 



(S(x, t)5 {n) (x - Xl )5 {m) (x-t-xl), (p{x, t 
= (H*A n H*B m ,S(x,t)<p(x,t)) 
= (H*(A n B m ),S{x,t)cp(x,t)) 
= - (A n B m , (Sip)(x + x u x - t + x ± - z*)) 
= - (5 {n \x) <g> 5 {rn) {t), (Sip)(x + x 1 ,x-t + x 1 -x* 1 ) 

= (-l) n+m+l d™d?(S V )(x + x u x - t + Xl - x\ 



x=0,t=0 
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n n+m 

= EE F 3i( x > t)&x&M x - t + X! - xl,X - t + Xl - Xl) 

j=0 i=0 
n n+m 

= EE F Jt (o,o)didM Xl ,ti) 

j=0 i=0 
n n+m 

= E E *ii(°,0) (^(a:-^)®^ ^-*!),^^) 

j=0 i=0 

Here Fji(x, t) are known smooth functions of 5 and of all its derivatives up to the 
order n + m. Hence, for all ip(t) G P(0, T) we get 

/ /" <5 (n) (a; - t)5 {rn \x -t-x\) dx^(t)\ 

n n+m I \ 

= E E **(°, °) ( / 6U) ( x - x i) ® <* w (* - *i) ^(*) ) 

j=0 i=0 \£ / 

n n+m 

= E E F ^(0, 0) - *J) ® ^(z - xx), l(x) ® m 

j=0 i=0 
n+m 

= £ f 0i (o,o) 

i=0 

As a consequence, 

I (t) — J b r (x)(S(x, t)a r (x — t) + Si(x, t)) dx 
t 

+(-l) n (S(x,t)a r (x-t) + S 1 (x,t))^\ x=xl (18) 

n+m 

+(-l) m d r x n S(x + t + x* 1 ,t)b r (x + t + x* 1 ) + V F 0l (0,0)5 {l \t-t* 1 ). 

Observe that the first three summands in (18) are smooth for t > 0. Indeed, 
the second summand is smooth due to a^(0) = for < i < n (see Assump- 
tion 1). The third summand is smooth due to b^\L) = for < % < m (see 
Assumption 2). 

Further plan of the solution construction. We split fif into subdomains 

n(i) = {(x, t) e ttj 1 1 - 1* < x < t - 

and construct the solution separately in each VL{%) and in a neighborhood of each 
border between Vt{i) and Vt{i + 1). Here £p = 0, 1 < i < k(T), where k(T) 
is defined by inequalities t* k ^ < T and t* k ^ +1 > T. The finiteness of k(T) is 
obvious. 
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6 Existence of the smooth solution on 



Lemma 21 There exists a smooth solution to the problem (8)-(10) on Q(l). 

Proof. Under the assumption that x\ < xi, we have t* < L. Hence (x\,tl) e 
Qq- Therefore any solution which is given by (13) on Qq, is smooth on fi(l), and 
has the property given by Item 9 of Definition 10, satisfies the integral Volterra 
equation of the second kind 

t—x 

u(x,t) = S 3 (x,t) + S 2 (x,t) J b r (£)u(Z,t - x) d£, (19) 



where 

S 2 (x, t) = S(x, t)c r (t — x) 

and 

S 3 (x,t) = S 2 (x,t)I (t — x) + S 1 (x,t) 

are known by (18). The smoothness of I Q (t — x) if (x,t) G follows from the 
facts that t — x < t\ and that I (t) restricted to the interval (0,t*) is smooth. 
Therefore S 2 and S% are smooth. 

The lemma will follow from two claims. Set 

Q*M(1) = {( x ,t) e £1(1) | i < t(m)}. 

Claim 1. Given m 6 No, there exists a unique solution u e C m (il*( m )(l)) to 
the problem (8) -(10) for some t(m) > 0. We apply the contraction principle to 
(19). Comparing the difference of two continuous functions u and u satisfying 
(19), we have 

\u — u\ < t(0)q max \u — u\, 

(a;,t)en*( )(i) 

where 

q— m ax \S\ max \c r \ max \b r \. 
(x,t)en(i) te[o,q] xe[o,L] 

Choosing t(0) < 1/q, we obtain the contraction property for the operator defined 

by the right-hand side of (19). The claim for m = follows. 

Our next concern is the existence and uniqueness of a C 1 (r2'( 1 )(l))-solution 
for some t(l). Let us consider the problem 

t—x 

d x u{x, t) = d x S 3 {x, t) + d x S 2 (x, t) J & r (£M£, 6(x, t)) d£ 



t—x 

-b r (t-x)u(t-x,t-x) - S 2 (x,t) J b r (£)(d t u)(£,t-x)d£. (20) 

o 
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From (8) we have d t u = p(x, t)u+g(x, t)—d x u. We choose an arbitrary t(l) < t(0). 
Since u is a known C(f2^ 1 - ) (l))-function, (20) on ^(^(l) is the Volterra integral 
equation of the second kind with respect to d x u. Assuming in addition to the 
condition t(l) < t(0) that t(l) < q, we obtain the contraction property for (20). 
On the account of (8), the claim for m — 1 follows. 

Proceeding further by induction and using in parallel (8), (19), and their 
suitable differentiations, we complete the proof of the claim. 

Claim 2. In the domain Q*i(l) there exists a unique smooth solution to 
the problem (8) -(10). Given m G N , we prove that there exists a unique 
u G C m (fi*i(l)) in at most \t\/t(m)\ steps by iterating the local existence and 
uniqueness result in domains 

n*'( m )(l)\n(*-i)'M(l), 1 < k < \T/t(m)]. 
In particular, for m = in the k-th step of the proof we have 

t-x-(k-l)t(m) 

u(x,t) = S 3 (x,t) + S 2 (x,t) J b r (0u(£,t - x) di 

^ 

+S 2 (x,t) J b r (0u(Z,t-x)d£ (21) 

t-x-(k-l)t(m) 

on {(x,t) G n kt(m \l) \x < t- (k- l)t(m)} 

and 

u(x,t) = S(x,t)u(0,t-x) + S 1 (x,t) on {(a;,*) G \ x> t- (k- l)t(m)}. (22) 

As in the latter formula t — x < (k — l)t(m), the function u defined by (22) is 
smooth and known from the previous steps. This implies that the last summand 
in (21) is known and smooth. Hence (21) is the Volterra integral equation of 
the second kind. Applying now the argument used to prove Claim 1, we obtain 
the existence and uniqueness of a continuous solution u to (21) on VL kt ( m ">(\) \ 
Q( fe_1 )*( m )(l). Since k is an arbitrary integer in the range 1 < k < \T/t(m)~\, we 
have u G C (Q**(l)). Further we similarly proceed with all derivatives f u. Claim 
2 is therewith proved. 

The solution on the whole Q(l) is now uniquely determined by the formula 

u(x, t) = S(x, t)u(0, t — x) + Si(x, t), 

where u(0, t — x) is a known smooth function. The latter is true due to < t — x < 
t\ and Claim 2. 

The proof of the claim is complete. □ 
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From the formulas (13) and (19), Lemma 21, and Assumption 1 it follows that 
u is smooth in a neighborhood of the characteristic line x = t. This ensures that 
u we construct satisfies Item 7 of Definition 10. 

Under the assumption that Q(2) is nonempty, in the next section we give the 
formula of the solution on 

n £ (l) = U {(x,t) G 0(2) \x > t-t\ -e} 

for a fixed e > such that i* — e > and 

6 r (x) = 0, xe[0,2e]. (23) 

Such e exists by Assumption 2. 



7 The solution on 

Write now 

L 

u(t) = J(b r (x)+5 (n \x-x 1 ))udx = v r (t)+v s (t), (24) 
o 

where v r (t) and t> s (t) are, respectively, regular (smooth) and singular parts of 
v(t). On the account of (16), (18), (23), and the fact that x\ < x±, we have on 
[0, *! + £]: 

t L 

v r (t) — J b r (x)u(x,t) dx + J b r (x) [S(x,t)a r (x — t) + Si(x,t)] dx 
t-q+e t 
+(-l) n d^(S(x,t)a r (x-t) + S 1 (x,t))\ x=Xl 
+ (-l) m 9™ (S(x + t + x{, t)b r (x + t + xl)) (25) 



x=0 



and 



n+m 
i=0 



(t)= ^F 0l (0,0)5«(t-t*)- (26) 



Note that the first summand in (25) is a known smooth function. This follows 
from the inclusion [t — £* + e,t] x {t} C 0(1) U {(x,t) \ x = t}, Lemma 21 and 
Assumption 1. 

We distinguish two cases. 

Case 1. t* = ti. As easily seen from (24), (25), and (26), v(t) = v r (t) on 
[0, tl+e] . Thus, Item 6 of Definition 10 for u we construct is fulfilled. Furthermore, 

u(0,t) = (5^(t-tl) + c r (t))v r (t) =J2CiS (i) (t-t*i) + c r (t)v r (t) (27) 

i=0 
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for t e (0, tl + e). The constants Cj depend on v^ k \t\) for < fc < j and can be 
computed by means of (25). 

Case 2. t\ ^ t 1 . Then x\ — x\ = t\. Using (24) and (25), we derive a similar 
formula for u(0, t) on (0, t\ + e): 

n+m n+m 

u(0,t) = c r (t) £ F 0i (0,0)5^(t-tl) + c r (t)v r (t) = E t 5^(t-tl) + c r (t)v r (t), 

i=0 i=0 

(28) 

where Ei are constants depending on F 0; fc(0, 0) and c^(ij) for < k < n + m. 
Set 



Q(0 = E c '< <j(i) (*-*D if^ = *i 

i=0 

and 

n+m 

Q{t)=Y J E l 5 {i \t-t\) iftl^h. 

i=0 

Lemma 22 u(x,t) given by the formula 

u(x, t) = S(x, t)c r (t — x)v r (t — x) + Si(x, t) + S(x, t)Q(t — x), (29) 

where v r (t) is determined by (25), is a V (Vl)- solution to the problem (8) -(10) 
restricted to fl e (l). 

Proof. On the account of (27), (28), and the construction of the solution on 
fi(l) done in Section 6, it is enough to prove that the restriction of S(x, t)Q(t — x) 
to Y — {0} x (0, t\ + e) is well defined and that S(x, i)Q{t — x) satisfies (8) with 
g(x,t) = on fi £ (l) in a distributional sense. The proof of the latter uses the 
argument as in the proof of Lemma 18. To prove the former, consider the smooth 
bijective map 

$ : (x, i) — > (x, t - x - t{). 

and its inverse 

: (x,t) -> (x,x + t + tl). 

Applying Theorem 6, we have 

WF($*flOc {(0,t + *;,-77, 77), 77^0}. 

Furthermore, 

N(Y) = {(p,t,t,0)} 

and therefore 

WF($* J B i ) n N(Y) =0 for all < i < n + m. 

By Theorem 4, the restriction of S(x, t)Q(6(x, t)) to Y is well defined. The lemma 
is therewith proved. □ 



21 



8 Construction of the smooth solution on Q(2) 



To shorten notation, without loss of generality we assume that maxf2f C\{(x, t)\x — 
0}>t* 2 . 

Lemma 23 There exists a smooth solution to the problem (8) -(10) on Q(2). 

Proof. We start from the general formula of a smooth solution on Q(2): 

u(x, t) = S(x, t)u(0, t-x) + S^x, t). (30) 

Since S and S\ are smooth, our task is to prove that there exists a smooth function 
identically equal to u(0, t — x) on fl{2). Since t\ < t — x < t 2 if (x, t) G Q(2) and 
c(t) = c r (t) if t G (tl,t 2 ), ^ suffices to show the existence of a smooth function 
iv(t) identically equal to v(t) on {t\,t* 2 ). From the formula (26) for v s (t) on 
(0, ^ + e) it follows that v(t) = v r (t) if t G (£[, + e), where e is as in Section 7 
and f r (t) is known and determined by (25). To prove the lemma, it is sufficient 
to show that there exists a smooth extension of v r (t) from (0, t\ + e) to [t* + e, t 2 ) 
such that v r (t) = v(t) if t G +e, i^)- If a suc h extension exists, then by (29) it 
satisfies the following integral equation on [£* + e,t* 2 ): 

t-q-s 

v r (t) = J b r (x)S(x, t)c r (t — x)v r (t — x) dx + R(t), (31) 



where 

p(t) P(t) 

R(t) — J b r (x)S(x,t)c r (t — x)v r (t — x) dx + J b r (x)Si(x, t) dx 

t-tl-e ( 32 ) 

+io(*) + J b r (x)S(x,t)Q(t - x) dx, 



6 r (x) is defined to be outside [0, L], and t> r in the formula (32) is known and 
defined by (25). One can easily see that the first three summands in (32) are 
smooth functions on [t*+e,t 2 ). We now show that the last summand is a C°°\t\ + 
e, ^-function as well. Indeed, take tp(t) G V(t\ + e/2,t\) and compute 



f b r (x)S(x, t)5 U) (t-x- tl) dx, ^(t) 
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= (5 {j \t - x - tl),b r (x)S(x,t)ij(t) 
= -(5 {j) (x) <g> l(t),b r (t - x - tl)S(t - x - tl,t)ijj(t)) 

= (-iy +1 (di (b r (t - x -ti)s(t-x- tit)) \ x=o , m 

We conclude that, irrespective of whether ti = t\ or t x ^ i*, the last summand 
in (32) is a known smooth function. As follows from (23), the functions v r {t) 
defined by (25) and (31) coincide at t = t\ + e. The same is true with respect to 
all the derivatives of v r . 

Our task is therefore reduced to show that there exists a C°°[t\+e, t^-function 
v r (t) satisfying (31). This follows from the fact that (31) is the integral Volterra 
equation of the second kind with respect to v r (t) (for details see the proof of 
Lemma 21). The proof is complete. □ 



9 Completion of the construction 

Continuing our construction in this fashion, we extend u over a neighborhood of 
each subsequent border between Vt(i — 1) and Vt{i) and over Q(i) for all 3 < i < 
k(T). Eventually we construct u on f2 T for any T > in the sense of Definition 10 
with f2 replaced by Vl T and II replaced by II T = {(x, t) G II 1 1 < T}. As easily 
seen from our construction, the condition (11) is fulfilled with Q + and Q' + replaced 
by Q T H Q + and fl T (1 respectively. Since T is arbitrary, the proof of Item 1 of 
Theorem 12 is complete. On the account of Definition 11 and the definition of the 
restriction u E V'(Q) to a subset of f2 (see [7, Section 5]), Item 2 of Theorem 12 
is a straightforward consequence of Item 1. Theorem 12 is therewith proved. 

Assume that S(x, t) ^ for all (x, t) e II. By (29) it follows from the 
construction, that if the singular part of b(x) is the derivative of the Dirac measure 
of order n, then for each % > 1 there exist j > i and n' > 1 such that u is the 
derivative of the Dirac measure of order n! along the characteristic line t — t* and 
u is the derivative of the Dirac measure of order n' + n along the characteristic line 
t — t*. In contrast, this is not so if singular parts of the initial and the boundary 
data are Dirac measures. In the latter case the solution preserves the same order 
of regularity in time. Furthermore, the assumption b^(L) = for all i £ N can 
be weakened to b r (L) = 0. 

Since u restricted to II \ I + is smooth, Theorem 16 follows from Item 2 of 
Theorem 15. 
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10 Uniqueness of the solution (Proof of Theo- 
rem 15) 

The proof of Theorem 15 is based on 5 lemmas. 

Lemma 24 A V' + {VL) -solution u to the problem (8) -(10) is unique on fi . 
Proof. Let u and u be two £>+(f2 )-solutions to the problem (8)-(9). Then 
(L(u-u),ip) = (u-u,L*ip) = for all ^6P(Q ), (33) 

where 

L = d t + \d x - p, L* = -(dt + \d x + p). (34) 
Our goal is to show that 

(u - u, -0) = for all ip G V(Q ). (35) 

Using the definition of T>' + (Qq) and (33), it is sufficient to prove that for every 
ip G V{VLq) there exists f G T>(Q ) such that 

L*(p = i/> on {(x,t) G Q \t > 0}. (36) 

Fix if; G P(fio)- If supp-0 fl {(x, t) 1 1 > 0} = 0, (35) follows immediately from 
the definition of T>' + (Q ). We therefore assume that supp-0 fl {(x, t) \ t > 0} ^ 0, 
Consider the problem 

<Pt + fx = - ^ (^ e {(x, t) g Q 1 1 > 0}, 

V^U=o = <Po(x), xe(0,L), 

where (fo(x) G V(0, L) will be specified below. This problem has a unique smooth 
solution given by the formula 

ip(x,t) = S(x,t)tp (x — t) + Si(x,t), 

where Si is given by (14) with p and g replaced by — p and —ip, respectively. 

Fix T(ip) > so that supp^ n {(x,t) \ t > T(ip)} = and S(x,T(ip)) ^ for 
all x with (x,T (■?/>)) G flo- The latter is ensured by Assumption 5. Set 

(/9 (X — -T('0)) = x 

S(x,T(^)) 
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for x such that (x,T(ip)) G IV Changing coordinates x — > £ = rr — T(^), we 
obtain 

m = gi(e + rW,rW) f37) 

5(e + T(^),T(^))" 
We construct the desired function <p(x,t) by the formula 

fq if (x,t) g g fi 1* > T($)}, 

ip(x, t) = < S(x, t)(p (x -t) + Si(x, t) if (x, t) G {(x, t) G Q | < t < T(ip)}, 
[(fi(x,t) if (x,t) G {(x,t) G tt \t < 0}, 

where (p(x,t) is chosen so that ip G V(Q, ). The proof is complete. □ 
From now on we use a modified definition of 

Q(i) = {(x, t) G n | * - t* < x < t - i > 1. 

Recall that t* = 0. 

Lemma 25 A V + (Q) -solution to the problem (8) -(10) is unique on 

Proof. Assume that there exist two £>+(f2)-solutions u and u. We will show 
that 

(v(t) - v(t),ip(t)) = for all ip(t) G P(0, £[), (38) 

where i> (t) is defined by Item 5 of Definition 10 and v(t) is defined similarly with 
u replaced by u. Postponing the proof, assume that (38) is true. Taking into 
account Item 2 of Definition 17 and the fact that c(t) — c r (t) if < t < t{, we 
have 

(L(u-u),(p) = {u-u,L*<p) = for all (p G V(Q(1)). 
Let us prove that 

(u-u,^) = for all tl)GD(fl(l)). (39) 

Following the argument used in the proof of Lemma 24, it is sufficient to show 
that, given ip G X>(fi(l)), there exists ip G X>(fi(l)) such that 

LV = ^ on {(x, t) G | x > 0}. 

We concentrate on the case that supp-?/> fl {(x,t) G f2(l) | x > 0} ^ 0. Otherwise 
(39) is immediate because u — w G X>+(fi(l)). Consider the problem 

V?t + ¥x = ~P<P - ip, (x, t) G {(x, t) G 0,(1) | x > 0}, 
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¥>U=o = <pi(t), t e (0,tj), 

where (fi(t) E X>(0, £*) is a fixed function. Let T(^) > be the same as in the 
proof of Lemma 24. We specify by 

= _A(T W -t.rW) (40) 
S(rw.)-{,rw)) 

and construct the desired ip similarly to the construction of tp in the proof of 
Lemma 24. To finish the proof of the lemma, it remains to show that 

(v -v,ip(t)) = for all ^(t) E V(ei, ei + 2e), (41) 

for each < i < t\/e — 2, where e > is chosen so that t\je is an integer and 

b r (x) = for x E [0, 2e]. (42) 

Such e exists by Assumption 2. We prove (41) by induction on i. 

Claim 1 (the base case). (41) is true for i = 0. We will use the following 
representations for u and u on Q + which are possible owing to Item 3 of Defini- 
tion 13: 

u = u Q + Mi in V(Q + ), , . 

u = u + iii in 

where uo = u and £to = w in X>'(Jlo H u = u = on Q 1 n Q + , U\ — u and 
Hi — u in D'^! fl Mi = "Si = on f2 fl Q + . 
We first prove that 

(v-v,rl;(t)) = {u 1 -u 1 ,b r (x)rl;(t)) for all #)eD(0,4e). (44) 

Accordingly to Item 1 of Definition 13, 

( V - V, ^(t)) = (U- Ub(x), 1(X) ® f45 A 

= ((u -u )6(a;),l(x)®^(t)) + ((ui-iti)6(a;),l(x)(8)V(t)), 1 ' 

where b r (x) — 0, x £ [0, L]. By Lemma 24, u = u in V(Qq fl Applying in 
addition Item 1 of Theorem 12 and Proposition 14, we have 

((«o - «o)(a;, t)b(x), l(x) ® ^(0) = (*(>(*) " *<>(*), ^(*)) , (46) 

where J (t) is defined by (17) and I (t) is defined by (17) with u replaced by u . 
From (18) we have Io(t) = Io(t) for < t < As. Hence the right-hand side of 
(46) is equal to 0. On the account of the inclusions supp(wi — ui) C VL\ and 
supp-^(t) C [0, 2e], (45) does not depend on b(x) outside [0,2s]. Since x* < xi, 
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b(x) = b r (x) on [0,2s]. Therefore (45) implies (44). Claim 1 now follows from 
(42). 

Assume that (41) is true ior i — k — 1, k > 1 and prove that it is true for 
i — k. 

Claim 2 (the induction step). (41) is true for i — k, k > 1. The proof is 
similar to the proof of Claim 1. Based on the induction assumption and applying 
the argument used in the proof of (39), we obtain 

u = u mV' + (G(k-l)), (47) 

where 

G(k) = £1(1) D{(x,t)\x>t-ek- 2e)}. 

Applying in addition Item 1 of Theorem 12, Proposition 14, and Lemma 21, we 
conclude that u is smooth on G(k — 1) n Q+. Owing to (47) and the latter fact, 
the following representations for u and u on Vt + are possible: 

u = u + Uk-i + Uk in V'(Q + ), 

u = u + u k _i + u k in V'(Q + ), 

where uq is the same as in (43), Uh-i = u in V(G(k — 1) n u k -i = on 
Q + \ G(k — 1), Uk = u and u k = u in V(Q + \ (G(k — 1) U O )) 7 u k = u k = on 
Vl + fl (G(k — 1) U Q ). Similarly to (44), we derive the equality 

(v — v, ip(t)) — (u k — u k , b r (x)ip(t)) for all ip(t) e V{ek, ek + 2e). 

The claim follows from the support properties of u k — u k , ip(t), and b r given 
by (42). 

The proof is complete. □ 
Set 

tt £ (0, 1) = {(x, t) e Q, | x - e < t < x + e)}. 

Lemma 26 A V + {VL) -solution to the problem (8) -(10) is unique on f2 £ (0, 1) pro- 
vided e is small enough. 

Proof. Let u and u be two V' + (r2)-solutions to the problem (8)-(10). Fix 
e > so that the condition (42) is fulfilled. By Claim 1 in the proof of Lemma 25, 
(41) is true for i — 0. Therefore 

(L(u-u),(p) = {u-u,L*(p) = for all ^GD(fi £ (0,l)). 
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Our task is to prove (39) with replaced by f2 e (0, 1). In fact, we prove that, 
given ip G V(£l e (0, 1)), there exists (p G V(il £ (0, 1)) satisfying the initial boundary 
problem 

ft + Vx = -p<p - ij>, (x,t) e fi e (o, l) n 

<p\t=o = <Po(x), x G [0,e), 
v\x=o = <pi(t), t e [0,5). 

Here tfo(x) £ C°°[0, e) is a fixed function identically equal to in a neighborhood 
of e, <fi(t) G C°°[0, e) is a fixed function identically equal to in a neighborhood 
of e, and ifo\o) = <Pi\o) for all i G No- We construct ip(x,t), combining the 
constructions of tp(x, t) in the proofs of Lemmas 24 and 25. Thus we fix T(ip) > 
to be the same as in the proof of Lemma 24 and specify ipo(x) and ipi(t) by (37) 
and (40), respectively. Let 

<p(x,t) 

f if (x, t) G {(x, t) G tr_(0, > T(^)}, 

t)v?oOr - t) + 5i(x, t) if (x, i) G {(x, t) G H^n fi £ (0, 1) | < t < T(ip)}, 
S(x,t)(pi(t-x)+Si(x,t) if (x,t) G {(x,t) G fi(l) nfi e (0, 1) |0 < t < T(i)}, 
{ (p{x,t) if (x,t) G {(x,t) G fi e (0, 1) \x < Oort < 0}, 

where <p(x,t) is chosen so that <p G £>(fi e (0, 1)). 

The proof is complete. □ 

For every % > 1 fix £j such that t* — £, > t* + £j < and 

b r (x) = for a; G [0,4^]. (48) 

Set 

Q(i) = {(;r, i) 1 1 - t* - Si < x < t - t* + Si}. 

Lemma 27 A V + (Q) -solution to the problem (8) -(10) is unique on Q(l). 

Proof. Assume that there exist two ^(f^-solutions u and u and show that 

(v-v,ip(t)) = for all ^(t) G V(t\ - £i,i* + £i). (49) 

By Lemmas 21 and 25, Item 1 of Theorem 12, and Proposition 14, any solution to 
(8)-(10) restricted to fi(l) is smooth. Based on this fact and on Lemmas 24-26, 
similarly to (44), we derive the equality 

(v - v, ip(t)) = (ui - «i, b r (x)ip(t)) 

for aRip(t) G T>(t\ — £±, t\ + £i), 
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where ui — u and u\ — u in V(G), u\ and u\ are identically equal to zero on 
n + \ G. Here 

G = {(X, t) E tt + | X < t - tl + £■!)}). 

The equality (49) now follows from the support properties of U\ — Hi, ip, and b r 
given by (48) for % — 1. 

We further distinguish two cases. 

Case 1. t\ ^ t x . Then c(t) = c r {t) for t in the range t\ — e 1 < t < t{ + E\. 
Applying (49) and Item 2 of Definition 13, we have 

L{u-u) = in V\Q{1)). (50) 

Case 2. t\ = h. Then c(t) = 5 {j) (t - ti) + c r (t). By Item 6 of Definition 10, 
v — v is smooth in a neighborhood of t\. Combining the latter with (49), we get 
(50). 

In the rest of the proof we proceed as in the proof of Lemma 25. □ 

Lemma 28 A V + (Q) -solution to the problem (8) -(10) is unique on f2(2). 

Proof. We follow the proof of Lemma 25 with replaced by f2(2) and 
with minor changes caused by the fact that due to Lemmas 27 and 23, u and u 
are smooth on Q(2) n f2 + H {(x, t) \x > t-t\-ei}. Hence (38) is true with (0, t\) 
replaced by (t\ + e 1 /2, t*). □ 

Continuing in this fashion, we eventually prove the uniqueness over subsequent 
Q(i) and Q(i) for any desired ? G N. Summarizing it with Lemmas 24 and 26 and 
Theorem 8, we obtain Item 1 of Theorem 15. 

Item 2 of Theorem 15 is a straightforward consequence of Item 1 of Theo- 
rem 15, Item 2 of Theorem 12, and Proposition 14. 
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